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ABSTRACT 

We present mult i- waveband optical imaging data obtained from observa- 
tions of the Subaru/XMM-Newton Deep Survey (SXDS). The survey field, cen- 
tered at R.A. = 02''18™00^ decl. = -05°00'00", has been the focus of a wide 
range of multi-wavelength observing programs spanning from X-ray to radio 
wavelengths. A large part of the optical imaging observations are carried out 
with Suprime-Cam on Subaru Telescope at Mauna Kea in the course of Sub- 
aru Telescope "Observatory Projects". This paper describes our optical obser- 
vations, data reduction and analysis procedures employed, and the character- 
istics of the data products. A total area of 1.22 deg^ is covered in five con- 
tiguous sub-fields, each of which corresponds to a single Suprime-Cam field of 
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view (~ 34' X 27' ), in five broad-band filters B,V, Rc,i' , z' to the depths of 
B = 28.4,1/ = 27.8, i?e = 27.7, z' = 27.7 and z' = 26.6 {AB,3a,(p = 2"). The 
data are reduced and compiled into five multi-waveband photometric catalogs, 
separately for each Suprime-Cam pointing. The «'-band catalogs contain about 
900,000 objects, making the SXDS catalogs one of the largest multi-waveband 
catalogs in corresponding depth and area coverage. The SXDS catalogs can be 
used for an extensive range of astronomical applications such as the number den- 
sity of the Galactic halo stars to the large scale structures at the distant universe. 
The number counts of galaxies are derived and compared with those of existing 
deep extragalactic surveys. The optical data, the source catalogs, and configu- 
ration files used to create the catalogs are publicly available via the SXDS web 



page ( jhttp : //www . nao j . org/Science/SubaruPro j ect/SXDS/index . htinl| ) . 



Subject headings: cosmology: observations — large-scale structure of universe 
galaxies: evolution — galaxies: formation — galaxies: photometry 



1. INTRODUCTION 

Understanding the formation and evolution processes of the individual galaxy and the 
growth of the large scale structures (LSSs) in the universe is one of the major goals in 
extragalactic astronomy today. In the scheme of a typical A-CDM model, the growth of 
structures is governed by the gravitational growth of initial fluctuations of dark matter. 
The baryonic material cools in these dark matter structures and grows through hierarchical 
clustering to galaxies and clusters of galaxies. The subsequent evolution of galaxies is closely 
connected to their environments, which, of course, relate to the LSSs, where those galaxies are 
located. The ultimate goal of extragalactic surveys is to trace these evolutionary processes 
by well-defined statistical galaxy samples. 

Optical imaging is arguably the cornerstone of any extragalactic survey, since it provides 
identifications and positions of celestial objects for follow-up spectroscopy. In recent years, 
there have been a number of deep imaging survey projects which devote significant amounts 
of telescope time, such as Great Observatories Origins Deep Survey (GOODS: Giavalisco et 
al. 2004) and Cosmic Evolution Survey (COSMOS: HST treasury project: Feldmann et al. 
2006; Scoville et al. 2007). These surveys provide multi- waveband galaxy samples at faint 



^ Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa- 
tory of Japan. 
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magnitude. Then, the photometric redshift techniques (Purusawa et al. 2000; Bolzonella et 
al. 2000; Feldmann et al. 2006) are frequently used to pre-select candidates for spectroscopy. 

Although a great deal of information can be obtained from the optical data alone, the 
value of the data set grows significantly as data at other wavelengths from other facilities 
are added. We therefore elected to use a multi- wavelength approach for the Subaru/XMM- 
Newton Deep Survey (SXDS; Sekiguchi et al. 2007, in preparation) from the very start, 
ensuring that our chosen field would be accessible and suitable for observations at all wave- 
lengths. An equatorial field, centered at R.A. = 02^18'"00^ dccl. = -05°00'00" is chosen to 
tic up with the deep X-ray observations with XMM-Ncwton observatory (Ueda et al. 2007), 
and also due to the accessibility by all major observatories, both existing and planned. The 
major multi-wavelength observation programs on the SXDS field (hereafter SXDF) include 
deep radio imaging with the VLA (Simpson et al. 2006), sub-mm mapping with SCUBA 
(Coppin et al. 2006; Mortier et al., 2005), mid-infrared observations with Spitzer Space 
Telescope (Lonsdale et al. 2003), deep near-infrared imaging with the UKIRT/WFCAM 
(Foucaud et al. 2006; Dye ct al. 2006; Lawrence et al. 2007), and the X-ray observations 
with XMM-Newton observatory. Importantly, the survey field of an infrared ultra deep 
survey (UDS; Foucaud et al. 2006) covering 0.77 square degrees as part of the UKIDSS 
project (Lawrence et al. 2007) is centered on the SXDF. It is expected that our extensive 
multi- wavelength data set will provide photometric redshifts accurate to 0.1 over a 

wide range of redshift, as well as detailed spectral energy distributions for the vast majority 
of the objects in the field. 

The SXDS has been undertaken as a part of the Subaru Telescope "Observatory Projects" , 
in which a large amount of observing times are devoted to carry out intensive survey pro- 
grams by combining observing times rewarded to builders and observatory staff of the Subaru 
Telescope. Note that Subaru Deep Field (SDF; Kashikawa et al. 2004) is another observatory 
project, which targets on a field different from the SXDF. The data set of the SDF is slightly 
deeper than the SXDF, but concentrates on one-fifth narrower field than the SXDF. Thus, 
these surveys complement each other and they can be used for a wide variety of studies. 

In this paper, we describe the observation details and data reduction and analysis of 
our optical imaging survey with the prime-focus camera on Subaru Telescope, as well as 
the resultant data products. A detailed description of the survey strategies, scientific ob- 
jectives, and multi-wavelength survey plans is given in a companion paper (Sekiguchi et 
al., 2007 - Paper I, in preparation). In Section 2, we describe optical imaging observations 
by Subaru Telescope, and in Section 3 we explain the data reduction procedure employed. 
We present creation of multi- waveband photometric catalogs including object detection and 
aperture photometry in Section 4, and calibrations of photometry and astrometry in Sec- 
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tion 5. In Section 6, we discuss characteristics of the catalogs such as hmiting magnitude, 
detection completeness, and number counts of galaxies. The summary is given in Section 7. 
Throughout this paper, all the magnitudes are expressed in the AB system unless otherwise 
mentioned. 



2. OBSERVATIONS AND DATA 

The optical imaging observations of the SXDF are carried out using the prime-focus cam- 
era (Suprime-Cam; Miyazaki et al. 2002) on Subaru Telescope in the period from September 
2002 to September 2005. Suprime-Cam is a 5 x 2 mosaic of ten MIT/LL buttable 2kx4k 
CCD camera. It has a 34' x 27' field of view with a projected pixel scale of 0.202 arcsec. The 
layout of the pointings is arranged as a cross shape so that each of the North-South and the 
East- West directions has an extent of ~ 1.3 deg (Figured]). This corresponds to a field span 
of a transverse dimension of ~ 75 Mpc at 2; ~ 1 and ~ 145 Mpc at ^ ~ 3 in the comoving 
scale for {h,QM,^A) = (0.7,0.3,0.7). The coverage of the comoving volume is ~ 3.2 x 10^ 
Mpc^ for z = — 3. Tabled] lists field center coordinates of the five pointings. Hereafter, the 
five pointings are referred to as C (Center), (North), S (South), E (East), and W (West), 
with respect to their relative positions on the sky. 

Observations of the SXDF are performed using five broad-band filters, B, V, Rc, i', and 
z' to cover the entire wavelength range observable with Suprime-Cam. The transmission 
curves of the five filters employed in the SXDS, convolved with the CCD sensitivities, reflec- 
tivity of the primary mirror, transmissions of corrector lenses, and atmospheric transmission, 
are given in Figure [2] Also, given in Figure [2] are the redshifted spectral energy distributions 
(SEDs) of an early-type galaxy at z = 0.5 and 1.2 and a late-type galaxy at ^ = 0.5 and 
5.0 for references. With photometric data sampled in the 5 bands, we can investigate color 
properties of red passively evolving galaxies to z ~ 1 (see Kodama et al. 2004; Yamada et al. 
2005) and young star-forming galaxies at high redshift such as lyman break galaxies (LEG) 
(Ouchi et al. 2004, 2005). The observations used in this paper span observing runs over a 
period of 3 years. A total of 160 hours are allocated to the SXDF imaging observations as a 
part of Subaru Telescope "Observatory Projects" (see also Kashikawa et al. 2004) combined 
with observatory staff times. In addition to these times, we use 14 hours of the observing 
times offered by the Supernova Cosmology Project (Doi et al. 2003; Lidman et al. 2005), 
which add the total exposure time of the center and the west pointings in the i' band {i'-C 
and i'-W). Moreover, in the course of another survey program for Lya emitting galaxies 
(Ouchi et al. 2005, 2007), 3 hours are devoted to the l^-band imaging. In total, 133 hours 
are used for on-source exposures for the final data set described in this paper. The complete 
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log of the observations is given in Table [2l 

The coordinates of the five pointings are carefully chosen so that the resultant images 
uniformly cover the entire SXDF, though we analyze each of the pointings separately in 
the present study. Also, to eliminate bad pixels, gaps between individual CCD chips and 
cosmic ray events, circular dithering pattern is used for each pointing. This dithering pattern 
employs a circular motion with radii of 60 to 120 arcsec around the center position of each 
pointing, combined with slight offsets to the center positions of the circular motion. As an 
example, the dithering pattern for the Rc band is plotted for each pointing in Figure O 

3. DATA REDUCTION 

The SXDS images obtained are processed using an in-house pipeline software package 
(See Ouchi 2003). We use it in almost the same way that successfully used for the SDF 
(Kashikawa et al. 2004). This pipeline is based on a software package developed for Suprime- 
Cam data reduction (Yagi et al. 2002) and several IRAF tasks (geomap, geotran, etc.) are 
used in the process of a geometrical transformation and alignment of final stacked images. We 
use a standard mosaic-CCD data reduction procedure, which is briefly summarized below. 

First, we apply a bias correction to the raw data. The median count of the overscan 
region is computed for each line, which represents a typical bias level at that line, and this 
median count is subtracted from the counts in all the pixels in that line. The overscan regions 
are trimmed after the bias subtraction. We assign a flag number (-32768) to the saturated or 
detected bad pixels, and these pixels are ignored in the following reduction processes. Since 
pixels in the outside edges of each frame are likely to be affected by noises which can not be 
corrected by the following fiatfielding process, they are also masked with the flag number. 

Flat frames are constructed from the normalized object frames taken during the same 
observing run. We use more than 150 object frames to create fiat frames. If the number 
of object frames is less than 150 frames, object frames taken in different observing runs 
within 3 months or twilight fiat frames taken in the same run are combined to increase the 
signal-to-noise ratios (S/Ns) of the fiat frames for that period. The median value of the 
each pixel is used to create the fiat frames. Pixels which could be affected by bright objects 
and vignetting by the auto-guider probe are flagged out before processing the flat frames. 
Then, object frames are divided by the flat frames. The peak-to-peak ratios of background 
levels over the entire flat fielded frames are less than 2%. We do not expect any significant 
systematic errors by the fiatfielding that affect accuracies of the following analyses. 

After the fiatfielding process is done, the distortion correction is applied to each frame 
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using the 5-th order polynomial transformation formula derived from Miyazaki et al. (2002). 

The changes in shapes and sizes of objects in the frames due to the distortion are negligible 
compared with the size of PSFs. Each pixel is transformed so that the surface brightness of 
the pixel is conserved, since the relative flux in each pixel has been already corrected by the 
fiatfielding process. Positional displacement of objects along the perpendicular orientation 
due to the atmospheric effect is also corrected at the same time. 

The sky background is subtracted from each distortion-free frame, We divide an image 
into temporary meshes with 64 x 64 pixels, corresponding to 12.9 arcsec in one direction. 
The mode of counts in each mesh is adopted as the sky background levels after smoothed by 
a median filter with a 3 x 3 kernel. This global sky subtraction should not affect photometric 
accuracies of compact objects such as galactic stars and faint galaxies. However, caution 
must be taken when we apply these values to extended objects such as nearby galaxies with 
low surface brightness. 

The stacking of the images is performed for each of the pointings separately in the 
following way. First, we choose 100 to 200 unsaturated stellar objects with a peak flux of 
> 500ADU in each frame distributing evenly over the entire frame. Second, we determine 
a reference frame of all the reduced frames. For the other frames, relative positions and 
rotation {Ax, Ay, A6) and flux ratios to the reference frame arc determined based on a least 
square method using these stellar objects commonly detected in the frames. Third, for all 
the frames, the positional shifts and rotations are corrected. Here, we do not use frames 
which have the flux ratios of less than 30% of that of the reference frame. Finally, all the 
frames are co-added by calculating flux-weighted average values in each pixel with a 3-a 
clipping process. The clipping process effectively removes cosmic rays and bad or saturated 
pixels, and satellite trails etc., and should not affect the resulting total fluxes of objects 
because the difference in counts among frames with various seeing sizes is mostly within the 
3-a threshold. 

After the co-adding of each band images, all the five-band {B, V, Rc, i', z') images are reg- 
istered into identical positional coordinates as we intend to execute multi- waveband aperture 
photometry. The positional registration is performed by the geometrical transformation by 
3-rd order polynomials fitting using the positions of stellar objects common to the co-added 
stacked images. The r.m.s. deviation of residuals of transformed positions with respect to 
the reference positions in the Rc-hand images is approximately 0.1-0.15 arcsec. 

For each pointing, the geometrically matched images are smoothed with Gaussian filters 

iteratively so that the PSF sizes of all the images become the same as that in the worst image 
used in this study. The PSF sizes of the resultant images are 0.80, 0.84, 0.82, 0.82, and 0.82 
arcseconds in SXDS-C, SXDS-N, SXDS-S, SXDS-E, and SXDS-W, respectively 
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Finally, regions with low S/Ns or those affected by saturation or overffow of electrons 
are carefully checked and assigned as flagged regions which are indicated in photometric 
catalogs. The resultant PSF sizes and total usable areas of each pointing is summarized in 
Table E 

4. MULTI- WAVEBAND PHOTOMETRIC CATALOGS 

We construct multi-waveband photometric catalogs of the objects identified in the 
SXDF. 



4.1. Object Detection 

Object detection and photometry are performed using SExtractor version 2.3.2 (Bertin 
& Arnouts 1996). Before performing the object detection and photometry, the outskirts 
of the stacked images are trimmed in order to eliminate any failure such as inappropriate 
estimation of sky levels and memory overflow during extraction of objects by SExtractor due 
to a high noise level. After trimming the outskirts, each individual image has similar size of 
approximately 10, 800 x 8, 100 pixels. We use the detection criteria that at least five pixels 
which have values of > 2a of sky background r.m.s. noise must be connected in order to be 
included in our object list. Positions of detected objects are translated to sky coordinates 
based on the astrometric calibration which will be discussed in Section 15. 2[ 

4.2. Aperture Photometry 

After the object detection is completed in one of the five bands, which we call 'detection 
band', multi- waveband aperture photometry is carried out with SExtractor for the other four 
bands using the same apertures at the same positions as those in the detection band. We 
perform aperture photometry with fixed circular apertures (MAG_APER) of 2 and 3 arcseconds 
in diameter, and variable elliptical apertures (MAG_AUTO and MAG_BEST) utilizing the Kron's 
first moment scheme with apertures' radii of 2.5 rxron (Kron 1980). Measured fluxes of objects 
are converted into magnitude using the photometric zeropoints, which will be described in 
Section 15.11 The extracted magnitude is considered to be the asymptotic total magnitude 
of galaxies, which should sample at least 94% of flux of a galaxy (Bertin & Arnouts 1996). 
Other useful parameters for objects including shape, size, and stellarity and so on, are also 
extracted by SExtractor at the same time. All the object parameters which are listed in the 
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catalogs are summarized in Table HI 

In this procedure, B-, V-, Rc-, i'-, and z'-selected multi-waveband catalogs are created 
for each pointing separately, which means that each object has measurements in the five 
bands {B, V, Rc, i', z') in each catalog. Thus 25 catalogs (= 5 pointings times 5 detection 
bands) in total are created. 

4.3. Flagged Areas 

Although we have fiagged out bad pixels from the image frame in the reduction proce- 
dure, there still remain (1) areas strongly affected mainly by overfiow of electrons or envelopes 
of extended objects, and (2) areas weakly affected by large halos from very bright stars or 
noisy regions. In our catalogs, the objects which are located in these areas are fiagged 1 and 
2, respectively. Otherwise, the objects located in the clean area have a fiag value of 0. 

The basic parameters used in the object detection and aperture photometry are summa- 
rized in Table [5l The complete configuration files for SExtractor can be found at|http: //www.naoj .org/Sci 
Note that we adopt default parameters of SExtractor for the source extraction unless we do 
not find any problem with the output values. The parameters adopted are suitable for ex- 
traction of moderate to faint compact galaxies, and not heavily optimized to extraction of 
extended objects or very faint objects. Thus the reliability of the current catalogs for those 
objects should be regarded with some caution. 

5. CALIBRATIONS 

5.1. Photometric Calibration 

Photometric zeropoints in all the five bands are determined in the following way. First, 
we compare our observation with the earlier calibrated observations for the same overlapping 
field, which were made by the Suprime-Cam instrument development team (Ouchi et al. 
2001, 2004). The southern half of our center pointing (SXDS-C) is overlapped with the 
northern half of the image obtained by Ouchi et al. (2001). They determined the photometric 
zeropoints in the Suprime-Cam band system based on observations of photometric standard 
stars SA92 and SA95 for B, V and Rc bands and a spectrophotometric standard star SA95-42 
for i' and z' bands. 

We estimate the photometric zeropoints for our center pointing [SXDS-C) by comparing 
magnitudes (MAG_AUTO) of ~ 200 objects which are commonly detected in the two images and 
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have FWHMs of smaller than 1.2 arcsec with a magnitude range of 20.0 to 24.0. Then, we 
extend these zeropoints to the surrounding four pointings by using objects in the overlapping 
areas in the same manner as for the determination of the zeropoint for the SXDS-C. In our 
catalogs, we list the instrumental magnitude obtained in the Suprime-Cam band system. 
Hence, no transformation of magnitude considering color terms of band responses are con- 
ducted. The color terms between the Suprime-Cam band system and the standard system 
are basically small - for the reddest stars, of order 0.1 in the B, V, i', z' bands, and slightly 
large (~ 0.3) in the Rc band. 

To test internal consistency of our zeropoint calibrations, we examine colors of stellar 
objects in our photometric catalogs and those computed by convolving typical stellar SEDs 
covering the complete ranges of spectral types (Gunn & Stryker 1983) with Suprime-Cam 
response curves. We find that it is necessary to apply corrections to the zeropoints of 
AV = +0.03 and Az' = +0.03 for all the five pointings, and Ai' = +0.05 and +0.02 for the 
South and East pointings {i'-S and i'-E), respectively. After these corrections are applied, 
the colors of the stellar objects in our catalogs agree with those computed for the SEDs from 
Gunn & Stryker (1983) within 0.03 mag in all the five bands (Figure HI). The systematic 
errors adjusted here may be due to relatively small number of stellar objects used in the 
determination of the photometric zeropoints. 

To ensure the accuracy of our zeropoints, we conduct CCD photometric observations of 
a part of the SXDF in B,Rc,i' and z' bands with the UH88-inch telescope at Mauna Kea 
on October 9 and 10, 2002. We compare magnitudes of objects in B,Rc,i' and z' bands 
obtained by the UH88-inch observations with our catalog values. We find that our catalog 
values agree within 0.05 mag r.m.s. for B and Rc bands and within 0.1 mag r.m.s. for i' 
and z' bands. The relatively large difference seen in i' and z' bands are probably due to 
the difference in the band responses between the UH88-inch system and the Suprime-Cam 
system, and low S/Ns of shallow SXDS images used for comparison with the UH88-inch data. 
To confirm accuracy of the photometric zeropoints in the i' and z' bands in our catalogs, 
follow-up observations are made for the two bands with the 1.0-m telescope at the United 
States Naval Observatory (USNO). The observations at the USNO were performed on the 
photometric night of December 3, 2004 in the course of photometric calibration for a high-^r 
supernovae search (Yasuda et al. 2007, in preparation). We find that a systematic difference 
between zeropoints in our catalogs and those determined by the USNO observation is as 
small as 0.03 mag r.m.s. 

Thus we conclude that the uncertainties of calibrated photometric zeropoints of the 
SXDS images are 0.03 — 0.05 mag r.m.s. The adopted photometric zeropoints are listed in 
Table M 
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5.2. Astrometric Calibration 

Astrometric calibration is performed using ~200 stars with magnitude of Jvega < 16.5 
derived from 2- MASS point source catalogs (PSC; Cutri et al. 2003), in each of the five 
stacked i?c-band images. Since images in the other bands have been geometrically trans- 
formed, the positions of objects detected on the images should coincide with those found in 
the i?c-band images (see Section [3]). The world coordinates for the images are calculated 
for the i?c-band images and written in the FITS headers of all the stacked images with the 
IRAF tasks CCMAP and CCSETWCS. The r.m.s. uncertainties of the determined coordi- 
nates across the FOV with respect to the PSC positions are on the order of 0.2 arcsec in 
both RA and Dec in each pointing (Tabled]). The positional error of the PSC is reported as 
~ 0.07 arcsec. 

We examine possible systematic differences in the calculated world coordinates between 
different pointings by using objects in the overlapping areas of two pointings. Figure [5] and 
Figure E] show the differences in RA and Dec of point-like objects with FWHMs < 1.2 arcsec 
in the magnitude range R = 20.5 — 24.0 detected both in the SXDS-C image and other 
surrounding images. From the figures, the world coordinates assigned to each of the five 
pointing images are in good agreement with one another with accuracy of about 1 arcsec 
at the edges of the images. The differences of the world coordinates between two pointings 
studied here are thought to be for the worst cases, since they are derived based on only 
the objects in the overlapping areas which are located at the edge of each image. Slopes 
and offsets in the residuals are seen in each panel of the figures. These disagreements in 
the positions of objects between two overlapping images come from possible tilts, offsets 
and/or difference in geometrical scale of the world coordinate systems between the two 
images. Residuals of distortion in the images and the atmospheric effect, which have not 
been completely removed in the reduction process, may cause the errors in determination 
of the world coordinate systems. Nevertheless, the SXDS catalogs have a good enough 
positional accuracy to perform follow-up spectroscopy. 

The astrometrically and photometrically calibrated images of the SXDS in five bands ( 
-B, V, i?c, i' , z') described in this section are released via SXDS Web Page (http : / /www . nao j . org/Science/J 
The data set described in this paper is based on the data release 1 (DRl). A complete history 
of the data release is summarized in the web site. The raw image data are also available to 
the public via Subaru Telescope ARchives System (STARS: https://stars.naoj .org). 
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6. PROPERTIES OF THE CATALOGS 

6.1. Galactic Extinction 

All the magnitudes listed in the catalogs are not corrected for the Galactic extinction. 
The magnitude attenuation in each band estimated for the central position (02'^18'"00'*, — 05°00'00"; 
J2000) based on Schlegel et al. (1998), assuming an extinction curve with Ry = Ay / E{B — 
V) = 3.1 is as follows: Ab = 0.091, Ay = 0.070, Ar^ = 0.056, A,/ = 0.044, and A.^ = 0.031. 
Table [7] summaries the magnitude attenuation in each band for all the pointings. The Galac- 
tic extinction above must be taken into account when studying extragalactic objects. 

6.2. Limiting Magnitudes 

We examine 3-a limiting magnitudes measured for the 2"-diameter apertures in all 
five bands {mum column in Table [2]) in the following manner. First, counts in unit of 
ADU fallen into the 2"-diameter apertures are measured at approximately 10,000 positions 
randomly selected and spreading over the entire image. Next, a histogram of the counts 
is produced and only the faint side of the histogram is fit by the Gaussian function, as 
the bright-side tail is composed of not only the sky background but also photons from 
the objects. The a of the best-fit Gaussian is regarded as the 1-a sky fiuctuation of the 
image for the 2"-diameter apertures. We find that limiting magnitude in each band is 
B = 28.4,1/ = 27.8, Rc = 27.7, i' = 27.7 and z' = 26.6 {AB,3(j,(p = 2") in the deepest 
images of the five pointings. 

6.3. S/N Distribution Map 

To examine the homogeneity of S/Ns of the SXDS data across the entire field, we 
investigate the S/N distribution map in the i?c-band images. First, we calculate the sky 
fiuctuation per 2"-diameter aperture within a mesh with a size of 350 x 350 pixels at each 
position in the same manner as in estimation of the limiting magnitudes. The mesh is shifted 
with a step of 175 pixels to cover the entire field of view. Then, the fiuctuations for the 2"- 
diameter aperture are converted to S/Ns for objects with an R^ magnitude of 27.5 at each 
position. 

Figure [7] shows the S/N distribution maps for 27.5-mag objects thus obtained for each 
pointing. The contour lines superimposed on the figures represent S/N=3 positions. First, 
we see that the brightest sources on the images significantly affect the S/Ns around the 
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sources. These areas are flagged in the catalogs. Second, each image shows almost an 
axisymmetric pattern of the S/N distribution when we ignore the low S/N areas due to the 
brightest sources. This pattern reflects the vignetting of the prime focus. Third, the pattern 
is slightly distorted for each image, probably due to bright sources, vignetting by the auto- 
guider probe, difference in the quantum efficiency among CCDs, and so on. However, the 
difference in S/Ns between in such distorted areas and in clean areas with higher S/Ns is 
on the order of 20% at the maximum. This difference is negligible in most studies even for 
faint sources. Thus, we can securely utilize the data with quite homogeneous S/Ns across 
the field of view. 



6.4. Completeness and False Detection Rate 

Completeness of the object detection as a function of magnitude for each band and for 
each pointing is estimated. The detection completeness is determined based on a Monte 
Carlo simulation by adding artificial objects which have Gaussian profiles with their Poisson 
photon noise into a stacked image used to create the catalogs and then detect them again 
in the same manner as in creating catalogs. The detection rate thus obtained is regarded 
as the detection completeness. Here the artificial objects are added to random positions 
in the effective area of the image. The FWHMs of the Gaussian profiles are set to be the 
same as the PSF sizes representative of the images which are derived from Table [3l The 
artificial objects contaminated by any neighbor sources listed in the catalogs are removed in 
calculation of the detection rate. If such blended objects are not removed from the sample, 
the detection rate should slightly decrease. However, since the contamination by neighbor 
objects leads to misidentification of the artificial objects in the object detection process, we 
adopt the approach which excludes the blended artificial objects, which was also employed 
in Kashikawa et al. (2004). 

The completeness determined by this procedure is shown in Figure [8] for all of the five 
bands. For the five pointings, the curves of the detection completeness drop similarly with 
increasing magnitude. We can summarize that the detection completeness in each image is 
at approximately 50% for objects with the 5-cr magnitude and 30% for objects with the 3-a 
magnitude, with a slight difference depending on the band. 

On the other hand, rates of the false object detection are also estimated as a function of 
magnitude. We generate negative images of each pointing and in each band by multiplying 
all the counts of the stacked images by —1, and then perform object detection for the negative 
images in the same manner as in creating the catalogs. In this process, detected objects are 
considered to be spurious objects. Thus we define the false detection rate as the number of 
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the spurious objects divided by the total number of objects in the catalogs in each magnitude 
bin, i.e., Ai'spurious("^)/^totai("^) (Figure[9]). It is seen that a contribution of the false detection 
is negligible, which does not exceed 0.5% if any, in the magnitude range brighter than the 
3-0" limiting magnitude. 

6.5. Magnitude Differences among the Catalogs 

In the areas where two images overlap with each other, a large fraction of objects are 
detected and their magnitudes are measured in both of the pointings. We can check system- 
atic difference in magnitude of those objects between the two catalogs. Since photometric 
zeropoints of the catalogs are determined so that MAG_AUTO of objects in the overlapping 
areas should coincide between any two catalogs (see Section [5TT|) . the magnitude difference 
in MAG_AUTO is negligible. MAG_BEST should be identical to MAG_AUTO for isolated objects. 
Therefore magnitude differences only for the 2"- and 3"-diameter aperture magnitudes are 
investigated in the following manner. 

First, we choose only compact objects with FWHMs < 1.2 arcsec which are listed in 
both the z'-selected catalog for SXDS-C and that for another pointing which overlap with 
each other, then measure the difference in magnitude of the objects by subtracting magnitude 
in the SXDS-C catalog by that in the other. Then, we determine the offset in the magnitude 
difference, i.e., systematic magnitude difference, by fitting the magnitude difference for a 
range of 21.0 to 23.5 mag by a least square method. Table [8] lists the systematic magnitude 
difference thus obtained for the 2"- and 3"- diameter aperture magnitudes in each band for 
each pair of catalogs. We note that 2"- and 3"-diameter aperture magnitudes indicate small 
systematic differences of < 0.05 mag between the catalogs. We think that this difference is 
probably due to difference in the PSF shape of the stacked images among the pointings. 

In order to see the effects of the systematic difference in the aperture magnitudes, we plot 
differences in colors of objects between the catalogs measured with the 2"-diameter aperture 
magnitude (Figure [TUl) . In this figure, for each of the four pairs of catalogs, the differences 
in colors {B — V^V — Rc, Rc — i', i' — z') are plotted. We see the systematic difference in the 
2"-diameter aperture colors of objects, which is no larger than 0.1 mag, due to the difference 
in the aperture magnitude between the catalogs. Since the magnitudes in the catalogs are 
not adjusted for this difference, the magnitude difference between the catalogs should be 
taken into account in the case that the 2"- or 3"-diameter aperture magnitude is used. For 
MAG_AUTO and MAG_BEST, which are considered as asymptotic total magnitude, the magnitude 
difference is negligible. 
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In the same manner as mentioned above, the r.m.s. of the magnitude differences for 
each band are calculated as a function of MAG_AUTO, which shows no systematic difference. 
Here the detection band for catalogs are chosen to be the same as the band for which the 
magnitude difference is investigated. For instance, to obtain the magnitude difference in the 
B band, the S-selected catalogs are used. Again, only the objects with FWHMs < 1.2 arcsec 
are used for measurement of the magnitude difference. The r.m.s. of magnitude differences 
measured in the above way are shown in Figure [TT] for each band. It is seen that the r.m.s. 
steadily increases with increasing magnitude. The r.m.s. of the magnitude difference can be 
converted into random magnitude error by dividing it by ^/2. We find that the magnitude 
errors thus estimated are entirely consistent with those expected from our calculation of 
the limiting magnitude discussed in Section 16. 2[ Note that the magnitude errors listed in 
our catalogs are computed by SExtractor assuming the simple Poisson statistics for the sky 
background fluctuations. The sky background noises measured in photometry apertures for 
the stacked images are likely to be larger than those estimated by SExtractor. Therefore 
the magnitude errors in the catalogs are underestimates and are not suitable for immediate 
analysis such as SED fitting and photo- 2; estimation. 

6.6. Number Counts of Galaxies 

To examine the characteristics of the catalogs we compare the number counts of galaxies 
with data from previous surveys. 

We calculate the number counts of galaxies by correcting our raw counts of detected 
objects against the detection completeness estimated in Section 16. 4[ The detailed studies 
of density fluctuations of galaxies based on the number counts or luminosity functions of 
galaxies can be found in elsewhere (e.g., Yamada et al. 2005). Star/galaxy separation 
is conducted based on FWHM and 2"-diameter fixed aperture magnitude of the objects, 
for which the effect of the Galactic extinction (Table [7]) is corrected. The criteria for the 
separation is determined so that both objects which have sizes as same as PSF of the images 
and saturated objects are effectively removed. For fainter magnitudes (> 22.5 — 23 mag), no 
separation is performed since a large fraction of galaxies have small sizes similar to those of 
stars, and contribution to the number count by stars is small at the faint magnitudes. The 
criteria for the separation adopted is summarized in Table O 

The results of the number counts of galaxies are as follows: First, internal comparison 
of the number counts among the five pointings of the SXDS are shown for each band in 
Figure [12] to Figure [161 In each figure, raw number counts of the galaxy samples extracted 
by the star/galaxy separation process above are plotted with open symbols. Effective ar- 
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eas used to calculate the number counts per unit area are given in Table [31 The number 
counts of galaxies which are corrected for the detection completeness are also superimposed 
with filled symbols. Here, the correction is performed by using the detection complete- 
ness as a function of magnitude which has been estimated for the Gaussian profiles, i.e., 
-?V(m) corrected = N (m) observed / completeness{m) . Error bars are calculated on the assumption 
of simple Poisson noise. We can say that no large systematic difference is seen in the superpo- 
sition of the number counts of galaxies among the five pointings in B,V, and Rc bands down 
to the 3-0" magnitudes. Similarly, in the other two bands, the number counts are consistent 
with one another among the five pointings within the Poisson errors at magnitudes brighter 
than ~ 25 mag. However, differences in the number counts at fainter magnitudes are found 
by a factor of 1.4 in the i' band and 1.7 in the z' band. This might be partly due to the 
field-to-field variation in the number counts of galaxies. However, since we have used only 
the Gaussian profiles with FWHM~ 0.8 arcsec in the simulation of the detection complete- 
ness, the completeness is likely to be overestimated at a faint magnitude range where most 
objects have extended profiles. So, we cannot conclude that the differences in the number 
counts at the faint end seen in the i' and z' bands indicates the field-to-field variation. 

Next, the number counts in the SXDF are compared with results by previous surveys 
(Figure[T7|to Figure [2T|) . Here, we calculate mean number counts of galaxies corrected for the 
detection completeness for the entire SXDF using the above sample for the five pointings, and 
plot them with filled circles with error bars. The error bars are again based on the Poisson 
noise. The results in other surveys are derived from Kashikawa et al. (2004), Capak et al. 
(2007), and Metcalfe (2007)@ and the references cited therein. A complete list of references 
compared in the figures is as follows: SDF - Kashikawa et al. 2004, HDFs/WHDFs - 
Metcalfe et al. 2001, Arnouts99 - Arnouts et al. 1999, SDSS - Yasuda et al. 2001, CADIS - 
Huang et al. 2001, Kummel& Wagner 01 - Kiimmel & Wagner 2001, Tyson88 - Tyson 1988, 
ArnoutsOl - Arnouts et al. 2001, Capak04 - Capak et al. 2004, Smail95 - Small et al. 1995, 
SH93 - Steidel & Hamilton 1993, Hogg97 - Hogg et al. 1997, COSMOS - Capak et al. 2007, 
Postman98 - Postman et al. 1998, Lilly91 - Lilly, Cowie, & Gardner 1991. 

From the comparison of the number counts among various surveys, we understand that 
our sample is consistent with results in other surveys for blank fields for the following reasons. 
First, the mean number counts of the SXDS show a good agreement with those from previous 
surveys down to magnitude where the number counts start to turn off. Second, it reproduces 
the faint-end data points derived from Hubble Deep Field (HDF). Again, the detection 
completeness simulated with the Gaussian profiles might cause a slightly large uncertainty 



" http : //star-www . dur . ac . uk/'mn/pubhtml/counts/ counts . html 
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of the corrected number counts at the faint end in each band. Thus, we conclude that the 
SXDS catalogs in the present study have characteristics similar to those in previous surveys 
and can be applied to general statistical studies on faint objects ranging from Galactic objects 
to high-redshift galaxies. 

We would like to mention how the Poisson error on the number counts is decreased by 
increase of the survey area. Figure [22] shows a comparison of the Poisson error associated 
with the number counts of galaxies among three surveys which are based on different survey 
areas, namely the HDF- North (5.3 arcmin^), GOODS-North (170 arcmin^) and Suprime- 
Cam 1 FOV (918 arcmin^). The uncertainty in the number counts due to the Poisson error 
is shown with a pair of lines for each of the surveys. For reference, the number count of 
SXDS is superimposed with a thick solid line. From the figure, we clearly see that the 
Poisson error decreases with increasing survey area. In particular, at the bright magnitudes 
of 21 — 23 mag, a large uncertainty by a factor of > 3 in the number count is found for 
the HDF-North area (dotted lines). For the same magnitude range, the uncertainty for a 
field of view of Suprime-Cam (dot-and-dash lines) is 8 — 15 times as small as that for the 
HDF-North area and 1.5 — 2.5 times as small as that for the GOODS-North area. Even at 
the faintest magnitudes of > 27 mag, a significant difference in the uncertainty from one 
to another survey area is seen. Thus it is emphasized that the SXDS data, which consists 
of five pointings of Suprime-Cam field of view, is a useful data set for studies on celestial 
objects which dominate a large fraction of the number counts without suffering from the 
Poisson error and the field-to-field variation as well, e.g., a study of evolution of the LSSs. 

7. SUMMARY 

The optical imaging observations of the SXDS project are carried out using the Suprime- 
Cam on Subaru Telescope in B,V,Rc,i' and z' bands. The SXDF has a contiguous area 
coverage of ~1.2 deg^, which consists of five Suprime-Cam pointings. 

The photometric zeropoints are determined with an absolute accuracy of no larger than 
0.05 mag r.m.s. in the photometry. The r.m.s. of the astrometric accuracies across the field 
is on the order of 0.2 arcsec in both RA and Dec. The systematic differences in the adopted 
world coordinates between different pointings are within about 1 arcsec at the outer edge 
of the field of view. Thus, our SXDS catalogs have a good enough positional accuracies to 
perform follow-up spectroscopy. 

Multi-waveband photometric catalogs of detected objects are created for each band and 
for each pointing. Each of the catalogs contains more than a hundred sixty thousand of 



- 18 - 



objects. The catalogs have quite homogeneous S/Ns across the field. The achieved limiting 
magnitudes in each band are B = 28.4, F = 27.8, i?c = 27.7, z' = 27.7 and z' = 26.6 
{AB,3a,(f) = 2"). 

The detection completeness as a function of magnitude is estimated by a Monte-Carlo 
simulation assuming a Gaussian profile. The number counts of galaxies in the SXDF in 
each band are computed by correcting the detection completeness, which are consistent with 
each other among the five pointings with a slight difference at faint magnitudes. The mean 
number counts of galaxies averaged over the five pointings show a good agreement with 
results from previous surveys down to the faint-end magnitude. 

With the aid of the wide coverage area of SXDS data, the uncertainty of the number 
counts of galaxies due to the Poisson error is greatly decreased. It is emphasized that the 
SXDS data is extremely useful for pursuing studies on celestial objects spreading in a wide 
field without suffering from the Poisson error and the field-to-field variation. The SXDS 
catalogs can be applied to studies, scientific objectives of which range from the Galactic 
objects to the large scale structures of the universe. 

The optical data, the compiled photometric catalogs, and configuration files used to cre- 
ate the catalogs have been released to the public and can be retrieved from the public data 
archives server of National Astronomical Observatory of Japan, |http : //www . nao j . org/Science/SubaruPr 
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Fig. 1. — Three color {B-, R^- and z'-band) composite optical image of the Subaru/XMM- 
Newton Deep survey Field (SXDF). The area consists of five Suprime-Cam pointings, cor- 
responding to a 1.22 square degree field. 
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Fig. 2. — Band response curves of the SXDS Suprime-Cam imaging observations (upper 
panel) and the typical spectra of an early- type galaxy (labeled as "elliptical") and a late- 
type galaxy (labeled as "disk" ) at the different redshifts (lower panel) . 
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Fig. 3. — Dithering patterns for the Rc-hand observations for five pointings. The center 
(0, 0) position of each pointing corresponds the coordinates given in Table [TJ 
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Fig. 4. — Color-color plots for the stellar objects in the SXDF (red dots) and Gunn-Stryker 
stars (blue triangles). The SXDF colors agree with the colors of Gunn-Stryker stars to the 
accuracy of 0.03 mag or better. 
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Fig. 5. — Residuals in the RA (filled circles) and in the Dec (crosses) of the World-Coordinate 
values in the RA direction in 4 overlapping areas between the center pointing (SXDS-C) and 
the north {SXDS-N: upper left panel), SXDS-C and the south {SXDS-S: upper right panel), 
SXDS-C and the east {SXDS-E: lower left panel), and SXDS-C and the west {SXDS-W: 
lower right panel). Only the objects with the FWHM< 1.2 arcsec are plotted. 
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Fig. 6.— 



Same as Figure [5], but the residuals in the Dec direction. 




Fig. 7. — Contour maps of the sky S/N distribution in the Rc band. The deviation (in ADU) 
of sky background fluctuations per 2"-diameter aperture are measured and converted to the 
S/Ns for objects with a magnitude of 27.5 at each position. 
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Fig. 8. — Object detection completeness of the five pointings as a function of magnitude in 
each band. 
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Fig. 9. — False object detection rate for the five pointings as a function of magnitude in 
the five bands. The perpendicular dashed lines indicate the 3-sigma limiting magnitude 
measured in 2"-diameter aperture in the center pointing {SXDS-C). 
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Fig. 10. — Residuals in colors measured for the 2"-diameter aperture magnitudes of the 
objects common in two separate pointing catalogs. 
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Fig. 11. — Random magnitude difference as a function of magnitude in each band. The 
random magnitude differences are computed by comparing fiuxes of compact objects with 
FWHMs < 1.2 arcsec between the central pointing and the four surrounding pointings. 
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Fig. 12. — 5-band galaxy number counts for the five pointings. Open symbols are the raw 
numbers extracted using the star/galaxy separation process. Filled symbols are the numbers 
corrected for the detection completeness estimated with Gaussian artificial sources. 
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Fig. 13. — \^-band galaxy number counts for the five pointings. Symbols are the same as in 
Figure [H 



-34- 



10^ 



I 

tie 



10^ 



R,-E 



10^ ^ 



i 



10^ ^ I 



10^ 



20 



22 



i 



a « 9 i 

S * 



e" 

B- 



24 



26 



28 



R 



C AB 



Fig. 14. — Rc-hand galaxy number counts for the five pointings. Symbols are the same as 
in Figure [T^ 
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Fig. 15. — i'-band galaxy number counts for the five pointings. Symbols are the same as in 
Figure [H 
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Fig. 16. — 2;'-band galaxy number counts for the five pointings. Symbols are the same as in 
Figure [H 
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Fig. 17. — Comparison of the galaxy number count of the SXDS (blue filled circles) with 
other major optical imaging surveys. The SXDS values are corrected for the detection 
completeness and the mean value for the five pointings is plotted. 
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Fig. 18.— 



Same as Figure [171 but for V^-band data. 
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Fig. 19.— 



Same as Figure [171 but for _Rc-band data. 
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Fig. 20.— 



Same as Figure [171 but for i'-band data. 
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Fig. 21.— 



Same as Figure [171 but for 2;'-band data. 
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Fig. 22. — i?-band number count uncertainties of galaxies, due to the Poisson error, for the 
Hubble Deep Field (HDF), GOODS North field (GOODS-N), a single Suprime-Cam field 
(SPCAM) are plotted with the SXDS number count (sohd line). The uncertainty for each 
survey is represented by a pair of lines plotted in the figure. 
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Table 1. COORDINATES OF SXDS FIVE POINTINGS 



Pointing 


RA (J2000) 


Dec (J2000) 


Position Angle of FOV 


Center 


02^18"00".000 


-05°00'00."00 


North is up 


North 


02'*18"00".000 


-04°35'00."00 


North is up 


South 


02^18"00".000 


-05°25'00."00 


North is up 


East 


02'^19"*47^070 


-05°00'00."00 


East is up 


West 


02'^16"*12^930 


-05°00'00."00 


West is up 



Note. — Field names for each pointing and their positional coordinates. PA, the orienta- 
tion of the FOV on the sky are chosen so that the overhead time to rotate the position angle 
are minimized and effective area with high signal-to-ratios are is maximized. 



Table 2. SUMMARY OF OBSERVATIONS AND DATA OF THE SXDS 



Exp. Time 


PSF Size 


Zeropts for Images 








Band-Pointing (min) 


(arcsec) 


(magAsADU-i) 


(magAs) 




Date of Observations 



B-C 


345 


0.80 


34.723 


28.09 


197,317 


2002, 


Sep. 


29/30, Oct. 01 














2003, 


Nov. 


17 


B-N 


330 


0.84 


34.701 


28.39 


176,372 


2002, 


Nov. 


02/05/27, Dec. 01/07 














2003, 


Nov. 


17 


B-S 


330 


0.82 


34.706 


28.33 


189,916 


2002, 


Nov. 


02/04/05/27, Dec. 01/07 














2003, 


Nov. 


17 


B-E 


330 


0.82 


34.698 


28.06 


179,478 


2002, 


Nov. 


04/05, Dec. 01/07 














2003, 


Nov. 


17 


B-W 


330 


0.78 


34.716 


28.21 


197,770 


2002, 


Nov. 


04/05/27, Dec. 01/07 














2003, 


Nov. 


17 


V-C 


319 


0.72 


33.639 


27.78 


213,851 


2003, 


Oct. 


02/24/26 














2004, 


Jan. 


16, Oct. 09, Dec. 10/14/1 














2005, 


Jan. 


05 


V-N 


313 


0.80 


33.648 


27.65 


203,299 


2003, 


Oct. 


02/26, Nov. 17 














2004, 


Jan. 


16, Dec. 10/14/15 














2005, 


Jan. 


05, Sep. 28 


V-S 


321 


0.82 


33.643 


27.75 


186,545 


2003, 


Oct. 


02/24/26 














2004, 


Jan. 


16, Dec. 10/14/15 














2005, 


Jan. 


05, Sep. 28 


V-E 


291 


0.76 


33.639 


27.77 


192,951 


2003, 


Oct. 


02/26 














2004, 


Jan. 


16, Dec. 10/14/15 














2005, 


Jan. 


05, Sep. 28 


V-W 


293 


0.72 


33.649 


27.72 


205,915 


2003, 


Oct. 


02/26, Nov. 17 














2004, 


Jan. 


16, Dec. 10/14/15 














2005, 


Jan. 


05, Sep. 28 


Rc-C 


248 


0.76 


34.315 


27.57 


188,079 


2002, 


Sep. 


30, Oct. 01/07 














2003, 


Nov. 


17 


Rc-N 


232 


0.78 


34.276 


27.74 


174,365 


2002, 


Nov. 


01/09/30, Dec. 06 














2003, 


Nov. 


17 


Rc-S 


232 


0.74 


34.219 


27.67 


178,416 


2002, 


Nov. 


01/09/30, Dec. 06 














2003, 


Nov. 


17 


Rc-E 


232 


0.76 


34.259 


27.51 


173,586 


2002, 


Nov. 


02/09/30, Dec. 06 



Table 2 — Continued 



Exp. Time 


PSF Size 


Zeropts for Images 








Band-Pointing (min) 


(arcsec) 


(magABADU-i) 


(magAB) 




Date of Observations 

















2003, 


Nov. 


17 


Rc 


,-W 


232 


0.82 


34.247 


27.53 


186,648 


2002, 


Nov. 


02/09/30, Dec. 06 
















2003, 


Nov. 


17 


i' 


-C 


647 


0.78 


34.055 


27.62 


181,352 


2002, 


Sep. 


29/30, Nov. 01/02/05/09/27/29, Dec. 06/07 
















2003, 


Oct. 


20/21 


i' 


-N 


440 


0.76 


34.042 


27.66 


176,394 


2002, 


Sep. 


29/30, Nov. 01/02/09/29, Dec. 07 
















2003, 


Sep. 


22, Oct. 02/21 


i' 


-S 


309 


0.68 


34.046 


27.47 


187,791 


2002, 


Sep. 


29/30, Nov. 01/02/09/29 
















2003, 


Sep. 


22, Oct. 02 


i' 


-E 


368 


0.74 


33.986 


27.49 


175,404 


2002, 


Sep. 


29/30, Nov. 01/02/09/29, Dec. 07 
















2003, 


Sep. 


22, Oct. 02/21 


%'■ 


-W 


598 


0.82 


34.087 


27.58 


178,543 


2002, 


Sep. 


29/30, Nov. 01/02/05/09/27/29, Dec. 06/07 
















2003, 


Oct. 


20/21 


z' 


-C 


217 


0.70 


33.076 


26.57 


183,324 


2002, 


Sep. 


29/30, Oct. 01, Nov. 04/05 


z' 


-N 


252 


0.74 


32.278 


26.64 


163,324 


2002, 


Nov. 


04/05/10 
















2003, 


Sep. 


21/22, Oct. 01/02 


z' 


■-S 


184 


0.76 


32.258 


26.39 


167,779 


2002, 


Nov. 


04/05/10 
















2003, 


Sep. 


21/22, Oct. 01 


z' 


-E 


267 


0.74 


32.743 


26.49 


160,415 


2002, 


Nov. 


04/05/10 
















2003, 


Sep. 


21/22, Oct. 01/02 


z' 


-W 


311 


0.74 


32.776 


26.60 


167,748 


2002, 


Nov. 


04/05/10 
















2003, 


Sep. 


21/22, Oct. 01/02 



Note. — Explanations of columns: (column 1) Filter pass-band {B, V, Rc, i' , or z') and the pointed field (center = 0, north 
= N, south =S, east = E, and west =W). (2) Total exposure time in minutes. (3) The PSF size of the stacked image. (4) 
Photometric zcropoint. (5) A 3-sigma limiting magnitude in AB magnitude, measured by random sampling with 2"-diameter 
aperture. (6) Number of objects detected in the image. (7) Date of observations. 
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Table 3. IMAGE SIZE AND THE AREA COVERED BY THE SXDS STACKED 

IMAGES 



Pointing 


PSF Size 
(arcsec) 


Area Covered 
(arcmin^) 


EfTective Area 
(arcmin^) 


Center (C) 


0.80 


1041.9 


968.4 


North (N) 


0.84 


1011.7 


951.8 


South (S) 


0.82 


1008.4 


967.3 


East (E) 


0.82 


986.4 


917.3 


West (W) 


0.82 


1017.7 


913.9 


Total 




4404.4 


4057.0 



Note. — Properties of stacked images of the SXDS. The table lists pointing field names, 
final PSF sizes, total areas covered by the stacked images, and effective usable areas after 
excluding the masked regions in each pointing. Five images in different bands in each 
pointing have the same PSF size. The total and effective areas covered by the 5 pointings 
taking account of overlapping areas between each pointing are 4404.4 arcmin^ and 4057.0 
arcmin^, respectively. 



Table 4. LIST OF PARAMETERS IN THE SXDS OPTICAL IMAGING OBJECT 

CATALOGS 



Parameter Name 


Unit 


Description 


ID 




sequential unique ID of the object in the catalog 


X 


pixels 


X coordinate of detection 


Y 


pixels 


Y coordinate of detection 


RA 


h:m:s 


right ascension of detection 


Dec 


d:m:s 


declination of detection 


KRON_RADIUS 


pixels 


Kron radius vk defined in Bertin & Arnouts 1996 


AJMAGE 


pixels 


2nd order moment along the major axis 


BJMAGE 


pixels 


2nd order moment along the minor axis 


ELLIPTICITY 




isophotal weighted eUipticity of the object (=A_IMAGE/B_IMAGE) 


THETAJMAGE 


degree 


position angle of the major axis, counter-clockwised, 0.0=X axis 


ISOAREAJMAGE 


pixels 


area of the lowest isophote 


Isop hot alM ag ( B ) 


mag 


isophotal magnitude in B 


IsophotalMag(V) 


mag 


isophotal magnitude in V 


IsophotalMag(R) 


mag 


isophotal magnitude in Rc 


IsophotalMag(i) 


mag 


isophotal magnitude in i' 


IsophotalMag(z) 


mag 


isophotal magnitude in z' 


IsophotCorMag(B) 


mag 


corrected isophotal magnitude in B defined by Bertin & Arnouts 1996 


IsophotCorMag(V) 


mag 


corrected isophotal magnitude in V defined by Bertin & Arnouts 1996 


IsophotCorMag(R) 


mag 


corrected isophotal magnitude in Rc defined by Bertin & Arnouts 1996 


IsophotCorMag(i) 


mag 


corrected isophotal magnitude in i' defined by Bertin & Arnouts 1996 


IsophotCorMag(z) 


mag 


corrected isophotal magnitude in z' defined by Bertin & Arnouts 1996 


2.0ApertureMag(B) 


mag 


2" fixed-aperture magnitude in B 


2.0ApertureMag(V) 


mag 


2" fixed-aperture magnitude in V 


2.0ApertureMag(R) 


mag 


(p 2" fixed-aperture magnitude in Rc 


2.0ApertureMag(i) 


mag 


(f) 2" fixed-aperture magnitude in i' 


2.0ApertureMag(z) 


mag 


(j) 2" fixed-aperture magnitude in z' 


3.0ApertureMag(B) 


mag 


3" fixed-aperture magnitude in B 


3.0ApertureMag(V) 


mag 


3" fixed-aperture magnitude in V 


3.0ApertureMag(R) 


mag 


3" fixed-aperture magnitude in Rc 


3.0ApertureMag(i) 


mag 


3" fixed-aperture magnitude in i' 


3. OApertureMag (z) 


mag 


3" fixed-aperture magnitude in z' 
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Table 4 — Continued 



Parameter Name 


Unit 


Description 




mag 


automatic-aperture magnitude in B 


A/TAP ATTTnA/"! 


mag 


duxomdxic-dperxure mdgnixuae in v 




mag 


automatic-aperture magnitude in Rc 


MAP ATTTn/'i'\ 


mag 


duxomdxic-dperxure mdgnixuae in i 


A/TAP \T\TD(r,\ 


mag 


duxonidxic-dpeixuic mdgnixuae in z 


A/TAP R'RCITf'R'l 


mag 


iVi/\vor_/\u i^j^^ijj 11 xiiere dre no neignuoio, or ibopnoxv_'Oiividg|^r>j oxne 


iVi/i.(ot_r>ll/o 1 V J 


mag 


MAG_AUTO(V) if there are no neighbors, or IsophotCorMag(V) othe 


iVl/\Vor_r>Il/0 i \iX) 


mag 


MAG_AUTO(Rc) if there are no neighbors, or IsophotCorMag(R) oth 


A/TAP R"I7QT(';'\ 


mag 


MAG_AUTO(i') if there are no neighbors, or IsophotCorMag(i) other\ 




mag 


iviAVj_f\u nj(^z j 11 xnere are no neignoors, or isopnoxL^oriviagi^zj oxner 


T7'yy^To/-\-rvVi/-\-i-ol A/To rf ^ /'R'\ 

Hirr ^^if5opnoT;dllVidg ) j 


mag 


isopnoxdi mdgnixuae rms error m d 


T^^■r^•(^arkr\\^rki-c^^\Ac^a^\ (\^\ 
HjL 1 ^^ibUpilULdllVidg )\y) 


mag 


ibopnoLdi nidgniLuae iiiib eiioi m v 


Hirr ^^if5opnoT;dllVidg ) ysx ) 


mag 


ibopnoxdi mdgnixuae rmb error m ric 


Err(IsophotalMag) (i) 


mag 


isophotal magnitude rms error in i' 


Err(IsopliotalMag) (z) 


mag 


isophotal magnitude rms error in z' 


Sill I l^ibOpilOLOOl iVidg J yD J 


Illdg 


coiiecxeQ isopiioxdi nidgiiixuGC iins eiioi in d 


H/I I l^inOpilOLOOl iVidgJ V j 


Illdg 


coiiecxeQ isopiioxai nidgiiixuQe iins eiioi in v 


Err(Isoph.otCorMag) (R) 


mag 


corrected isophotal magnitude rms error in Rc 


Hirr (^ibopiiotv^orividg jyi) 


mag 


correcxea isopnoxai mdgnixuae rms error in i 


Illrr^^if5opnoT;^^orlVidgJ yz j 


mag 


correcxea isopnoxdi mdgnixuae rms error in z 


Illrr^^z.u/\perxureiVidgJ^^IJ ) 


mag 


nxea-dperxure mdgnixuae rms error m jj 


ii/rr(^z.UApertu.reiViagj(^ V ) 


mag 


02" fixed-aperture magnitude rms error in V 


Err(2.0ApertureMag) (R) 


mag 


02" fixed-aperture magnitude rms error in Rc 


Err(2.0ApertureMag) (i) 


mag 


02" fixed-aperture magnitude rms error in i' 


Err(2.0ApertureMag) (z) 


mag 


02" fixed-aperture magnitude rms error in z' 


Err(3.0ApertureMag) (B) 


mag 


03" fixed-aperture magnitude rms error in B 


Err (S.OApertureMag) (V) 


mag 


03" fixed-aperture magnitude rms error in V 


Err(3.0ApertureMag) (R) 


mag 


03" fixed-aperture magnitude rms error in Rc 


Err(3.0ApertureMag) (i) 


mag 


03" fixed-aperture magnitude rms error in i' 


Err(3.0ApertureMag) (z) 


mag 


03" fixed-aperture magnitude rms error in z' 


Err(MAG_AUTO)(B) 


mag 


automatic-aperture magnitude rms error in B 
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Table 4 — Continued 



Parameter Name 



Unit Description 



Err(MAG_AUTO)(V) mag 

Err(MAG_AUTO)(R) mag 

Err(MAG_AUTO)(i) mag 

Err(MAG_AUTO)(z) mag 

Err(MAG_BEST)(B) mag 

Err(MAG_BEST)(V) mag 

Err(MAG_BEST)(R) mag 

Err(MAG_BEST)(i) mag 

Err(MAG_BEST)(z) mag 

FLUX_MAX(B) ADU 

FLUX_MAX(V) ADU 

FLUX_MAX(R) ADU 

FLUX_MAX(i) ADU 

FLUX_MAX(z) ADU 

FWHMJMAGE(B) pixels 

FWHMJMAGE(V) pixels 

FWHM_IMAGE(R) pixels 

FWHM_IMAGE(i) pixels 

FWHM_IMAGE(z) pixels 
CLASS_STAR(B) 
CLASS_STAR(V) 
CLASS_STAR(R) 
CLASS_STAR(i) 
CLASS_STAR(z) 
FLAGS(B) 
FLAGS (V) 
FLAGS(R) 

FLAGS(i) 
FLAGS(z) 
FLAG_AREA 



automatic-aperture magnitude rms error in V 
automatic-aperture magnitude rms error in Rc 
automatic-aperture magnitude rms error in i' 
automatic-aperture magnitude rms error in z' 
MAG_BEST magnitude rms error in B 
MAG_BEST magnitude rms error in V 
MAG_BEST magnitude rms error in Rc 
MAG_BEST magnitude rms error in i' 
MAG_BEST magnitude rms error in z' 
peak surface brightness in B 
peak surface brightness in V 
peak surface brightness in Rc 
peak surface brightness in i' 
peak surface brightness in z' 
FWHM of profile from a gaussian fit in B 
FWHM of profile from a gaussian fit in V 
FWHM of profile from a gaussian fit in Rc 
FWHM of profile from a gaussian fit in i' 
FWHM of profile from a gaussian fit in z' 
stellarity index in B: 0.0=Galaxy and 1.0=Star 
stellarity index in V: 0.0=Galaxy and 1.0=Star 
stellarity index in R: 0.0=Galaxy and 1.0=Star 
stellarity index in i': 0.0=Galaxy and 1.0=Star 
stellarity index in z': 0.0=Galaxy and 1.0=Star 
extraction fiags in B defined by S Extractor 
extraction flags in V defined by SExtractor 
extraction flags in Rc defined by SExtractor 
extraction fiags in i' defined by SExtractor 
extraction fiags in z' defined by SExtractor 

index for affected area: 0=Clean, 1=STR0NG_HAL0, 2=WEAK_HAL( 
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Note. — The list of fields listed in the optical multi-waveband catalogs. All the fields are 
commonly included in all the catalogs created for each detection and each pointing. 
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Table 5. THE SEXTRACTOR PARAMETERS USED TO CREATE THE SXDS 

CATALOGS 



Parameter 


Value 


DETECT_MINAREA 


5 


DETECT_THRESH 


2.0 


ANALYSIS_THRESH 


2.0 


FILTER 


N 


DEBLEND_NTHRESH 


32 


DEBLEND_MINCONT 


0.005 


BACK_SIZE 


32 


BACK_FILTERSIZE 


3 


BACKPHOTO_TYPE 


GLOBAL 


BACKPHOTO_THICK 


24 


SATURXEVEL 


50000.0 


GAIN 


2.6 



Note. — Parameters for SExtractor used to create catalogs. The above values are common 
to all the SExtractor runs regardless of detection bandpasses. 
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Table 6. RESULT OF ASTROMETRIC CALIBRATION 



Pointing 


rms (RA) 
(arcsec) 


rms (Dec) 
(arcsec) 


Number of Stars 


Center 


0.160 


0.151 


233 


North 


0.220 


0.218 


208 


South 


0.251 


0.207 


194 


East 


0.261 


0.225 


199 


West 


0.217 


0.201 


218 



Note. — The resultant r.m.s. of astrometric calibration. The number of stars used for 
calculation of the world coordinates is shown in the far right column. 



-53- 



Table 7. THE GALACTIC EXTINCTION IN THE SXDF 



Pointing 


E{B - V) 


Ab 


Av 


Arc 


Ai' 


A,, 


SXDS-C 


0.021 


0.091 


0.070 


0.056 


0.044 


0.031 


SXDS-N 


0.019 


0.083 


0.064 


0.052 


0.040 


0.029 


SXDS-S 


0.023 


0.100 


0.077 


0.062 


0.048 


0.034 


SXDS-E 


0.021 


0.091 


0.070 


0.056 


0.044 


0.031 


SXDS-W 


0.021 


0.091 


0.070 


0.057 


0.044 


0.031 



Note. — The magnitude attenuation due to the Galactic extinction in each band calculated 
for the central coordinates of each pointing based on Schlegel et al. (1998). Since magnitudes 
in the released catalogs are not corrected for the Galactic extinction, the above values must 
be taken into account when studying extragalactic objects. 
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Table 8. THE SYSTEMATIC DIFFERENCES IN APERTURE MAGNITUDES 
BETWEEN THE DIFFERENT CATALOGS OF THE SXDS 



Pointing 


B{2") 


S(3") 


V{2") 


y(3") 






i'{2") 


i'(3") 


z'{2") 


z'{?>") 


C-N 


-0.04 


-0.02 


-0.04 


-0.01 


-0.03 


-0.01 


+0.03 


+0.02 


-0.02 


-0.01 


C-S 


-0.05 


-0.02 


-0.02 


-0.00 


-0.03 


0.00 


+0.03 


+0.00 


+0.00 


+0.01 


C-E 


-0.02 


-0.01 


-0.02 


-0.01 


-0.01 


-0.00 


+0.05 


+0.01 


-0.01 


-0.00 


C-W 


+0.01 


+0.01 


-0.03 


-0.01 


+0.01 


+0.00 


-0.01 


-0.00 


-0.01 


-0.01 



Note. — The systematic differences in aperture magnitudes of objects between the cat- 
alog of SXDS-C and that of another surrounding pointing. For each band, the magnitude 
differences in 2"- and 3"-diameter aperture magnitude are listed. The column of pointing 
denotes the pair of pointings compared. For instance, the line for C-N indicates magnitude 
differences calculated by magnitudes in SXDS-C minus those in SXDS-N. 
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Table 9. CRITERIA USED FOR STAR/GALAXY SEPARATION 



Band 






Criteria 










B 


{m = 


19.5 


- 23.0, a > 4.5,(7 > -2.5 (m 


-20.3) + 4.5) 


or 


(m 


> 23.0) 


V 


{m = 


19.5 


- 23.0, cT > 4.5,(7 > -2.5 (m 


- 20.0) + 4.5) 


or 


(m 


> 23.0) 


Rc 


{m = 


19.5 


- 22.5, (T > 4.5,(7 > -2.5 (m 


-20.3) + 4.5) 


or 


(m 


> 22.5) 


i' 


{m — 


19.5 


- 22.5, (7 > 4.5, (7 > -2.5 (m 


-20.1) + 4.5) 


or 


(m 


> 22.5) 


z' 


(m = 


18.5 


- 22.5,(7 > 4.5,(7 > -3.0 (m 


- 19.3) + 4.5) 


or 


(m 


> 22.5) 



Note. — The criteria used in star/galaxy separation procedure. Objects satisfying these 
criteria are recognized as galaxies. In the criteria list, m denotes 2" -diameter fixed aperture 
magnitude and a is FWHM of objects in the corresponding band. 



